Intensity dependence of multiple orbital contributions and shape resonance in 
high-order harmonic generation of aligned N 2 molecules 



O 

(N 
•*-> 
O 

O 

00 



CZ5 

O 

q 
o 

>V 

Oh. 



> 

o 



X 



Cheng Jin, 1 Julien B. Bertrand, 2 R. R. Lucchese, 3 H. J. Worner, 4 
Paul B. Corkum, 2 D. M. Villeneuve, 2 Anh-Thu Le, 1 and C. D. Lin 1 

1 J. R. Macdonald Laboratory, Physics Department, 
Kansas State University, Manhattan, Kansas 66506-2604, USA 
2 Joint Laboratory for Attosecond Science, National Research Council of Canada and University of Ottawa, 
100 Sussex Drive, Ottawa, Ontario K1A 0R6, Canada 
s Department of Chemistry, Texas A&iM University, College Station, TX 77843-3255, USA 
* Laboratorium fur Physikalische Chemie, ETH Zurich, 
Wolfgang-Pauli-Strasse 10, 8093 Zurich, Switzerland 
(Dated: October 19, 2011) 

We report measurements and theoretical simulations of high-order harmonic generation (HHG) in 
aligned N2 molecules using a 1200-nm intense laser field when the generating pulse is perpendicular 
to the aligning one. With increasing laser intensity, the minimum in the HHG spectra first shifts its 
position and then disappears. Theoretical simulations including the macroscopic propagation effects 
in the medium reproduce these observations and the disappearance of the minimum is attributed 
to the additional contribution of HHG from inner orbitals. We also predict that the well-known 
shape resonance in the photoionization spectra of N2 should exist in the HHG spectra. It is most 
clearly seen when the generating laser is parallel to the aligning one, and disappears gradually as 
the angle between the two lasers increases. No clear evidence of this shape resonance has been 
reported so far when using lasers with different wavelengths. Further experimentation is needed to 
draw conclusions. 

PACS numbers: 33.80.Rv, 42.65.Ky, 31.70.Hq, 33. 80. Eh 



I. INTRODUCTION 



Recently high-order harmonic generation (HHG) has 
been employed to probe the electronic structure of 
molecules on an ultrafast time scale dri]. After an elec- 
tron is first removed from the molecule by the infrared 
(IR) laser, it may be driven back by the action of an in- 
tense laser field to recollide with the parent ion. HHG is a 
process where the returning electrons recombine with the 
molecular ion with the emission of high-energy photons. 
In fact, the last step - - "recombination" is the inverse of 
photoionization. Any spectral features in the photoion- 
ization cross section (PICS) would thus be embodied in 
the HHG spectra as well. Since molecules can be impul- 
sively aligned by a laser field @, 0], the observation of 
HHG from aligned molecules further offers the opportu- 
nity for probing PICSs from aligned molecules that are 
not generally possible with synchrotron radiation experi- 
ments. In addition, the first step of HHG is a highly non- 
linear tunneling ionization process which is very selective 
with respect to the ionization energy of the molecular or- 
bital from which the electron is removed. Thus HHG in 
general is dominated by the recombination to the highest 
occupied molecular orbital (HOMO). In contrast, direct 
photoionization of a molecule is a linear process, where 
PICSs from the next or next few inner orbitals could 
be of comparable importance. For example, for N2, the 
HOMO is a a g orbital. The next more tightly bound 
orbital, the HOMO-1, is a ir u orbital with an ionization 
energy 1.3 eV higher. The HOMO-2, is a a u orbital and 
is 3.2 eV more tightly bound than the HOMO. Calcu- 



lations by Lucchese et al. |8[ and experimental results 
0, [l(| show that for photon energy from the threshold 
at 15 eV to about 40 eV, the PICS from the HOMO-2 is 
negligible, but the HOMO-1 PICS is actually comparable 
with that of the HOMO. On the other hand, the HOMO 
has a shape resonance that peaks near the photon en- 
ergy of 29 eV. These are predictions and measurements 
made on randomly oriented N2 molecules. How do these 
features depend on the alignment of molecules? Today 
the rich structure in PICSs from fixed-in-space molecules 
predicted by the theory remains mostly unexplored ex- 
perimentally. Can high-harmonic spectra generated by 
laser pulses from aligned molecules provide new informa- 
tion that are not yet directly available from photoioniza- 
tion measurements? The answer is yes. Already experi- 
ments [ill [T2I ] have shown the existence of minimum in 
the HHG spectra from aligned molecules which has not 
been seen in photoionization experiment. Moreover, Mc- 
Farland et al. [13] reported that HHG from the HOMO-1 
dominates over that from the HOMO when the molecules 
are perpendicularly aligned with respect to the laser po- 
larization. Mairesse et al. [Tij performed harmonic spec- 
troscopy to characterize the attosecond dynamics of mul- 
tielectron rearrangement during strong-field ionization. 

Earlier attempts to extract structure information from 
HHG spectra were based on the two-center interference 
model 15] and the strong-field approximation (SEA) fl6l — 
[Lsj . Subsequently, a quantitative rescattering (QRS) the- 
ory [l9l - [2~fl | was developed which established that HHG 
spectra from an isolated molecule can be expressed as the 
product of a returning electron wave packet with a photo- 
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recombination cross section (PRCS). The QRS asserts 
that the PRCS is independent of the laser parameters, 
including the wavelength, intensity, and pulse duration. 
The latter affect the returning electron wave packet only. 
In the meanwhile, Smirnova et al. 0] studied HHG from 
aligned CO2 molecules and emphasized the importance of 
hole dynamics with including multiple orbitals. Within 
the QRS theory, multiple orbitals can be easily incorpo- 
rated into the theory and it was first used by Le et al. [22[ 
to explain the HHG data by McFarland et al. [l3[ . These 
theoretical studies were all based on single-molecule cal- 
culations. To compare with experimental measurements, 
however, the effect of the propagation of the driving laser 
field and the harmonics in the medium has to be consid- 
ered. Until recently, the propagation effect has only been 
considered for atomic targets, but not for molecules [23l — 
|26| . Even for atoms, the propagation effect is considered 
mostly only within the SFA. Since experimental HHG 
spectra depend critically on laser parameters, focusing 
conditions, gas-jet pressure, and how the spectra are col- 
lected, these parameters have to be available in order 
to simulate the experimental spectra. When the experi- 
mental conditions are well specified, we have shown that 
direct comparison with experimental HHG spectra is now 
possible, for example, for N 2 [HHi] and C0 2 HI- These 
comparisons are also becoming more interesting with the 
emergence of the use of mid-infrared lasers to cover a 
broader range of photon energies. 

In this combined theoretical and experimental paper, 
we revisit the HHG spectra of N 2 that are aligned per- 
pendicular to the laser polarization. This was previously 
studied by McFarland et al. [l3[ using 800-nm lasers 
near the cutoff region. Here we use 1200-nm lasers. We 
revisit the issue about the importance of the HOMO-1 
contribution raised in Refs. [IH, [22], Wf\ . We find that 
HOMO-1 does not contribute to the HHG spectra at low 
laser intensities, but becomes comparable to the HOMO 
at higher intensities. We also examine the shape reso- 
nance in the differential PICS, and show how they would 
appear in the HHG spectra of aligned N 2 . The paper 
is organized as follows. In Sec. II, information on the 
experimental setup is given. In Sec. Ill, we briefly sum- 
marize the theoretical methods used in the simulation. 
In Sec. IV, we show the comparison between experimen- 
tal and theoretical HHG spectra, the differential PICS of 
the HOMO and HOMO-1, and the shape resonance in 
the HHG spectra of the HOMO for aligned N 2 . A short 
summary in Sec. V concludes this paper. 



II. EXPERIMENTAL SETUP 

To perform the experiment, we use the output of 
a Ti:Sapphire multi-pass laser system (32 fs, 800 nm, 
50 Hz, 12 mJ per pulse), an optical setup to both split 
and recombine laser pulses, and a high harmonic cham- 
ber, described previously [HI, composed of a source 
chamber (pulsed valve, 250 /itm orifice) and an extreme 



ultraviolet (XUV) spectrometer. 

First, we impulsively align molecules [H, 0] with a 
stretched non-ionizing pump pulse (70 ± 5 fs, 800 nm, 
^a/jgn=5xl0 13 W/cm 2 ). Second, we probe them with 
a delayed intense high harmonic generation pulse (1200 
nm: 40 ± 5 fs of variable intensity). The 1200 nm 
pulses (~ 1 mJ) come from a high-energy optical para- 
metric amplifier (HE-TOPAS) pumped by 800 nm light 
(~ 8 mJ). In both arms of our Mach-Zehnder interferom- 
eter, the intensity is adjusted with neutral absorption- 
type density filters. Both linearly polarized pump and 
probe beams are focused (f=50 cm) ~ 1 mm downstream 
in the gas jet expansion and ~ 2 mm before the orifice 
to select only the so-called "short" trajectories [29j | . 

We record high-harmonic spectra at maximal align- 
ment at the half revival (~ 4.12 ps). Our experimen- 
tal conditions (Pbacfe=2 atm., T rot =30-40K, I & u gn = 
5 x 10 13 W/cm 2 , and T a / ig „=70±5 fs) suggest we reach 
a maximum degree of alignment of (cos 2 6*}=0.60±0.05 
based on calculations [3(| and supported by recent su- 
personic gas expansion studies in similar conditions [3l| . 
We vary the angle a between the alignment axis and 
the probing field by rotating the alignment beam polar- 
ization, both beams are linearly polarized. This paper 
concentrates on the case where a — 90°. 



III. THEORETICAL METHODS 

The details of the theoretical methods used in the sim- 
ulations have been described in Refs. 0, HH. Briefly, 
both the fundamental laser field and the XUV light are 
modified when they co-propagate through a macroscopic 
medium. For some circumstances, such as low laser inten- 
sity, low gas pressure, and short gas medium, the effects 
of dispersion, Kerr nonlinearity, and plasmas defocusing 
on the fundamental laser field can be neglected 0, H2] • 
Under these conditions the profile of the fundamental 
laser field in space (in the vacuum) can be expressed in 
an analytical form. If the IR laser is considered to be a 
Gaussian beam, its spatial and temporal dependence is 
given approximately in Ref. (32j. For the propagation 
of the harmonics in the medium, the dispersion and ab- 
sorption effects, which are generally anisotropic in case 
of aligned molecules, are not included when the pressure 
is low [24]]. We only include the induced dipoles for the 
harmonics generated by the IR laser. The harmonics 
emitted from the gas medium propagate in the vacuum 
after the gas jet, then they are detected by the spectrom- 
eter. A Hankel transformation [33| is applied to obtain 
the far-field harmonics using the harmonics at the exit of 
the gas medium. 

The three-dimensional Maxwell's wave equation for the 
XUV light in the gas medium is 0, [3i437| 

v^H, , v l d 2 El(r,z,t,a) d 2 P^(r,z,t,a) 
V El(r,z,t,a) - ^ = W> ^ •(!) 

Here E^(r, z, t, a) and Pj^Ar, z, t, a) are the parallel com- 
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ponents (with respect to the polarization direction of the 
generating laser) of the electric field of the XUV light 
and the nonlinear polarization caused by the IR laser, 
respectively, a is the pump-probe angle, i.e., the angle 
between the aligning (pump) and the generating (probe) 
laser polarizations. 

The nonlinear polarization term can be expressed as 



oc=90° 



Pr!i(r,z,t,a) = [n - n e (r, z,t,a)]D n ' ° (r,z,t,a), 



(2) 



where tiq — n c (r, z, t, a) gives the density of the remaining 
neutral molecules, and D^ ,tot (t, a) is the parallel compo- 
nent of the induced single-molecule dipole over a num- 
ber of active electrons (including the effects from outer- 
most and the inner orbitals). The QRS theory |19l - [2l| is 
applied to calculate the induced dipole D^' tot (t, a) from 
each orbital. Within the QRS, the laser- induced dipole 
moment D(uj,9) for a fixed- in-space molecule (the molec- 
ular axis has an angle 9 with respect to the laser polar- 
ization) is given by 



£>H(w,0) = N(6) l/2 W(u)S(u,6), 



(3) 



where N(0) is the alignment-dependent ionization proba- 
bility, W(u>) is the recombining electron wave packet, and 
d"(oj, 9) is the parallel component of the photorecombi- 
nation (PR) transition dipole (complex in general). 

For partially aligned molecules, we define an ionization 
weighted PR transition dipole moment by 



d avg (w,a) 



N(9) 1 ^ 2 S (w, 9)p(9, a) sin 9d9, (4) 



where p{9, a) is the alignment distribution in the probe- 
laser frame [38| . Such a definition can be made for each 
orbital from which the ionization occurs. The harmonic 
signal from each orbital, after propagation, can be ex- 
pressed as m m m 
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FIG. 1: (Color online) Comparison of experimental [red (dark 
gray) lines] and theoretical [green (light gray) lines] HHG 
spectra of aligned N2 in a 1200-nm generating laser where the 
pump-probe angle a=90° . Laser intensities in the simulations 
are indicated where Io = 10 14 W/cm 2 . Degree of alignment is 
(cos 2 #}=0.60 by the aligning laser. See text for additional 
laser parameters. Only a orbital is included in the simula- 
tions in (a) and (b), and both a and n orbitals are included 
in the simulation in (c). Arrows indicate the positions of min- 
ima. 



IV. RESULTS AND DISCUSSION 



^( W ,a) KW 4 |W>)|V v ^,a)| 2 , (5) 

where W(u) is the "macroscopic wave packet" (MWP). 
If more than one molecular orbital contributes to the 
HHG spectra, then the transition dipole will be added 
up coherently. Note that in Eq. (|3|) the N(6) for each 
molecular orbital is calculated using the MO-ADK the- 
ory. 

Based on the previous studies of N 2 molecules, only 
the HOMO and HOMO-I orbitals, which are a and ir 
orbitals, respectively, need to be included in the HHG 
calculation. We will use either the terms HOMO and 
HOMO-I, or a and 7r orbitals, interchangeably below. 



A. Macroscopic HHG spectra: Experiment vs 
Theory 

The HHG spectra for unaligned and aligned N2 and 
CO2 molecules using 800- nm and 1200-nm lasers have 
been reported recently [39]. Comparison with theoreti- 
cally simulated spectra including macroscopic propaga- 
tion has been carried out by Jin et al. [231 - 1251 ]. In 
this paper we present new measurements for aligned N2 
molecules at a time delay corresponding to the maximal 
alignment while the HHG-generating laser is perpendic- 
ular to the aligning one. The degree of alignment is esti- 
mated to be (cos 2 6>}=0.60±0.05. 
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Fig. Q] shows the HHG spectra of aligned N2 generated 
by a 1200-nm laser which is perpendicular to the align- 
ing laser. Experimentally, the laser duration (FWHM) 
is ~ 44 fs, the beam waist at the focus is ~ 40 /im, and 
the length of gas jet is ~ 1 mm, which is located 3 mm 
after the laser focus. A vertical slit with a width of 100 
/im is placed 24 cm after the gas jet. The intensities (in 
the center of the gas jet) used in the theory are adjusted 
to coincide with experimental HHG cutoff position. In 
Figs. [Ha)- HJc), the laser intensities in theory (experi- 
ment) are 0.75 (0.65), 0.9 (1.1), and 1.1 (1.3), in units 
of 10 14 W/cm 2 , respectively. We use the degree of align- 
ment (cos 2 #}=0.60 and keep other parameters close to 
experimental ones in the simulation. The HHG spectra 
from the experiment and theory are normalized at the 
cutoff. 

The main features in the spectra are the deep minima 
at 38.2 eV and 40.4 eV, at the two lower intensities in 
Figs. QJa) andQJb), respectively. The minimum disap- 
pears at the higher intensity in Fig. HJc). To simulate 
the spectra at the two low intensities, we include the a 
orbital only. The simulation reproduces not only the cor- 
rect shape of the spectra, but also the precise positions 
of the minima in the spectra in Figs. [IJa) andQJb). For 
the higher intensity in Figs. [He), when we only included 
the a orbital the theory could not reproduce the correct 
spectral shape. It also predicted a minimum in the spec- 
trum which was not seen in the experiment [see Fig.[2Ja)] . 
We then included both a and it orbitals. A very good 
agreement between theory and experiment (correct shape 
and no minimum in the spectrum) in Fig. [IJc) is then 
achieved. 



B. Single HOMO orbital contribution at low laser 
intensity 

We now take a careful examination of the spectral fea- 
tures in Figs. Ufa) andQJb). The deep minimum is re- 
lated to the a orbital. This minimum has been observed 
in many experiments, either in unaligned or aligned N2 
[39TJ42l | . The minimum shifts only slightly when the laser 
intensity is changed. This behavior is similar to the well- 
known Cooper minimum in Ar [23], IH, 44] . The same be- 
havior with laser intensity has been observed by Worner 
et al. (see their Fig. 1) and Farrell et al. [42| (see 
their Fig. 7) when a=0°. This can be understood as 
proposed in Ref. [23jj. When only one molecular orbital 
is contributing to the HHG spectra, the harmonic yield 
is given in Eq. ([5]). In Fig. [5] we show that there is a fast 
drop-off in the averaged PR transition dipole (shown is 
the square of the magnitude multiplied by w 2 ) of the a 
orbital around 38 eV [see Fig. Hfb)], which causes the 
pronounced minimum in the HHG spectra. For different 
laser intensities, the slope of the averaged PR transition 
dipole does not change, but the MWP (cj 2 |VF(u;)| 2 ) in- 
creases monotonically with photon energy. Since the spa- 
tial distribution of the IR laser intensity is different when 
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FIG. 2: (Color online) (a) Calculated macroscopic HHG spec- 
tra (envelope only) corresponding to Fig. [He) • Total (HOMO 
and HOMO-1 together) spectra and individual HOMO and 
HOMO-1 spectra, (b) Averaged photorecombination transi- 
tion dipoles (parallel component, the square of the magnitude 
multiplied by uj 2 ) of the HOMO and HOMO-1 corresponding 
to (a). Laser intensity is l.lxlO 14 W/cm 2 . Degree of align- 
ment is (cos 2 t?}=0.60 by the aligning laser. The pump-probe 
angle a=90°. 



the peak intensity at the focus is increased, the wave 
packet is modified differently as it propagates through 
the medium (see Fig. 3(d) in (23[). This makes the min- 
imum in the HHG spectra shift slightly with laser inten- 
sity. Using only one molecular orbital, our theoretical 
simulation can reproduce the position of the minimum 
close to the experimental one. 



C. Multiple orbital contributions (HOMO and 
HOMO-1) at higher laser intensity 

We next analyze the spectral features in Fig. [He). In 
Fig. [3a), the envelopes of the HHG spectra in Fig. [TJc) 
from the two individual molecular orbitals, and from the 
total one, are obtained after macroscopic propagation in 
the gas medium, are shown. Meanwhile, the averaged 
PR transition dipoles (the square of the magnitude mul- 
tiplied by uj 2 , degeneracy is not included) of the two or- 
bitals are shown in Fig.^b). 

In Fig.^a), the a orbital alone shows a deep minimum, 
similar to Figs. (Ha) andQJb). Since the tt orbital shows 
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comparable contributions over the whole spectral region, 
the interference between a and tt orbitals washes out the 
minimum in the spectra. From Eq. (j4|) , the relative con- 
tribution between the two orbitals can be adjusted by the 
alignment-dependent ionization probabilities N(9) and 
the alignment distribution p(9,a). For lower intensities, 
the relative N(6) of the tt orbital is small, so only the 
HOMO orbital contributes to the spectra. As the in- 
tensity increases, both orbitals contribute and the two 
amplitudes interfere resulting in a drastic change of the 
spectra. 

As shown in Eq. (O, the macroscopic HHG spectra 
from individual molecular orbitals can be considered as a 
product of a MWP and an averaged PR transition dipolc, 
based on the QRS theory. Since the ionization potential 
of the a orbital (15.6 eV) differs from the tt orbital (16.9 
eV) only by 1.3 eV, the MWPs of the two orbitals are 
almost the same under the same IR laser, so the rela- 
tive contribution between a and tt orbitals to the total 
HHG spectra is mostly determined by the averaged PR 
transition dipoles. In Figs. HJa) and [U(b), the magni- 
tude of N(9) for a orbital is much larger than the one for 
the tt orbital, thus making the corresponding averaged 
PR transition dipole also larger. At the higher inten- 
sity of Fig. Q7c), the averaged transition dipoles between 
the two orbitals become comparable. Thus by increas- 
ing the laser intensity, the total HHG can evolve from 
single-orbital to multi ple- orbital phenomena. Note that 
the MO-ADK theory [H, 51] is used to calculate N[&}. 
There are also other models in the literature [22, El 



for calculating the ionization rates. Using different ion- 
ization rates and different alignment distributions may 
change the theoretical predictions. We find that the ion- 
ization rates obtained from MO-ADK theory is very close 
to the recent model calculation used by Petretti et al. 
[48j where they solved the time-dependent Schrodinger 
equation at laser intensity of 1.5xl0 14 W/cm 2 . 

The contribution of the tt orbital to the HHG of N2 
when t he pump -probe angle a=90° has been studied pre- 
viously [13], [23] using 800- nm pulses at the high intensities 
of around 2.0x 10 14 W/cm 2 . In this case, HOMO-1 was 
found to be much more pronounced in the cutoff region. 
Note that the QRS theory has been applied to interpret 
the results in 13l| with the ionization probability obtained 
from SEA [22| without including the propagation effect. 
(Note that at the alignment angle of 90°, the ratio of 
ionization rate between HOMO and HOMO-1 is about 2 
from the SFA in comparison with about 5 from the MO- 
ADK theory (H| Ef3].) In this study, a 1200-nm laser 
generates a broad photon energy range even with a low 
laser intensity and the tt orbital contributes not only in 
the cutoff region, but also in the plateau. 

In the future, measurements similar to the present one 
but with a full range of pump-probe angles may provide 
a way to determine the relative ionization probabilities 
of the two orbitals. The pump-probe angle a=90° is 
much closer to the alignment angle of 90° since the "vol- 
ume element" sin 8d9 in the angular integration peaks at 



90°. To probe the PR transition dipole of the HOMO 
orbital, HHG spectra taken at low laser intensity with 
long-wavelength lasers is preferable to avoid multiple or- 
bital contributions. This would make the retrieval of the 
structure of the target easier. If one wishes to study the 
tt orbital, then a higher laser intensity and better align- 
ment will enhance its contribution to the HHG [13, . 



D. Photoionization amplitudes and phases of N2 
from HOMO and HOMO-1 orbitals 

The most basic information on photoionization of a 
fixed-in-space molecule is the parallel and perpendicu- 
lar transition dipole amplitudes in the body-frame of the 
molecule. They appear in the differential PICS in the 
body-fixed frame as [1, |49| 



4ir 2 ujk 



dQ;dQn 



1 f,k' 



(6) 



where h is the polarization direction of the light, k the 
momentum of the photoelectron, and uj the photon en- 
ergy. We will only focus on the case of h \\ k since they 
are related to the parallel polarized HHG spectra mea- 
sured from aligned molecules. 

In Fig. [3] we show the differential PICSs (in units of 
Mb) and the corresponding phases of the a and tt or- 
bitals of No using the well-established photoionization 
theory @, E<| for photon energy from 20 to 80 eV. These 
data have been shown previously for selected photon en- 
ergies against the alignment angles for the HOMO (see 
Fig. 5 of UJ) and HOMO-1 (see Fig. 2 of [22]). The ob- 
served HHG minima shown in Fig.[T]above for a=90° and 
in Fig. 2 of [23[ for a=0° and for randomly distributed 
N2 molecules can all be understood as due to the rapid 
change of cross section near 40 eV at alignment angles 
either close to 0° or 90°. Clearly the precise position 
of the minimum will depend on the degree of alignment. 
For the HOMO-1, Fig. [3J shows that the cross section 
generally peaks at large alignment angles. Thus interfer- 
ence of HHG spectra from the a and tt orbitals are only 
observed close to a=90°. Note that in Fig. [3] the PICS 
for the two orbitals are shown on the same scale. Except 
for the HOMO shape resonance near 30 eV that will be 
discussed next, the PICS from HOMO-1 is comparable 
with that from HOMO over the photon energies covered. 
This result has been confirmed by the measured electron 
spectra in photoionization experiments. In HHG, HOMO 
contribution is always dominant for the randomly dis- 
tributed N2. HOMO-1 can become more important only 
at large alignment angles and higher laser intensities, as 
discussed earlier. 



E. Shape resonance in HHG of aligned N2 

Resonances are very common in photoionization spec- 
tra. Most of them are due to the so-called Feshbach 
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FIG. 3: (Color online) Calculated differential photoionization cross sections [(a) and (c)] and phases [(b) and (d)] of the HOMO 
and HOMO-1 for N2 in terms of the alignment angles, respectively. The shape resonance in the HOMO shows up around 30 
eV for small alignment angles only. Only the parallel component to the polarization direction of the laser is shown. 



resonances which in general are quite narrow and thus 
can only be observed using high-resolution spectroscopy. 
However, broad Feshbach resonances and shape reso- 
nances can be observed since they have widths from frac- 
tions to a few eV's. In HHG, resonances have been ex- 
plored in experiment (Hoj and theory [5lH54| recently for 
atomic targets. But shape resonances are rare for com- 
mon atomic targets. The resonance feature observed in 
Xe is due to the intershell coupling with the well- 
known shape resonance that occurs in the 4d shell. On 
the other hand, shape resonances are very common in 
molecules. For N 2 , there is a pronounced shape resonance 
in the HOMO channel near photon energy of 30 eV. This 
is due to the 3<j g —tkcr u channel, and for small alignment 
angles only. We can see a decrease of the phase shift by tt 
for this resonance from 20 to 40 eV in Fig.[3][b). There are 
no other known shape resonances in the covered energy 
region. Clearly this shape resonance is best observed by 
selecting ionization from the a orbital only, by using low 
laser intensity and for molecules that are aligned nearly 
parallel to the polarization axis of the probing laser. 

In Fig. SI we first show the calculated HHG spectra 
(envelope only, normalized) of randomly distributed N2 
at two laser intensities. For randomly distributed N2 
molecules, it is known that the a orbital is the dominant 
contributor to the HHG spectra. The peak around 30 eV 
in the HHG spectra is due to the shape resonance in the 
PICS of a orbital at small alignment angles. In Fig. [H 
we also show the normalized yield at selected alignment 
angles with respect to the randomly distributed one. We 



take the degree of alignment (cos 2 6*)=0.60. The inten- 
sity of each odd harmonic q is obtained by integrating 
over harmonics q — 1 to q + 1. For low laser intensity, 
the shape resonance is very pronounced at a — 0°. It 
decreases as the pump-probe angle is increased, showing 
that the shape resonance is present only at small align- 
ing angles. For higher laser intensities, the same behavior 
is seen even though the tt orbital is also contributing at 
large pump-probe angles. We should comment that ab- 
sorption is not included in the present simulation and it 
will suppress the shape resonance if the gas pressure is 
high. 

Comparison of this prediction with existing experi- 
mental data is far from conclusive. Torres et al. (40| 
have shown high harmonics data for aligned N2 using 
a 1300-nm laser with intensity of 1.3x 10 14 W/cm 2 (see 
their Fig. 4) . The general trend of their data is very close 
to our prediction, but they used a higher intensity and a 
higher degree of alignment ((cos 2 6*}=0.66), and thus the 
HOMO-1 channel may contribute to the signal at larger 
alignment angles. Using 800-nm lasers, Lee et al. |27| 
reported the HHG ratios of aligned vs unaligned N2 at 
the selected alignment angles. They did not present data 
near the resonance region and they used a high intensity 
of 2.5xl0 14 W/cm 2 [see their Fig. 1(c)]. Our calcula- 
tion does not reproduce their measured ratios. They also 
reported that the normalized HHG yields for a=0° and 
90° crosses at about harmonic order 39. This is con- 
firmed in our calculation (not shown). The very recent 
measurement by Kato et al. (5o| did not extend below 
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FIG. 4: (Color online) Calculated macroscopic HHG spec- 
tra (envelope only) of unaligned N2 and the normalized HHG 
yields of aligned N2 with respect to the unaligned ones as a 
function of the photon energy at the selected pump-probe an- 
gle a. Laser intensities are indicated where Io=10 14 W/cm 2 . 
Degree of alignment is (cos 2 9) =0.60. 



may modify the prediction if the gas pressure is too high. 
Experiments dedicated to address this issue would be of 
interest. 



V. CONCLUSIONS AND OUTLOOK 

In this paper, we have reported experimental HHG 
spectra of aligned N2 at a time delay of the maximal 
alignment using a 1200-nm laser when the pump-probe 
angle is set at 90°. The minimum of the HHG spectra 
appears at about 38 - 40 eV at two low laser intensities. 
It disappears at a higher laser intensity. We have carried 
out theoretical simulations to understand these results, 
and concluded that the minima in the HHG are associ- 
ated with the properties of the transition dipole moments 
for photoionization from the HOMO. At higher intensity, 
the contribution from HOMO-1 becomes important and 
interference between the two contributions washes out 
the minimum. 

We have also examined the possibility of observing the 
well-known shape resonance in the photoionization of N2 
in the HHG spectra. While the normalized HHG yield 
(with respect to randomly distributed molecules) shows 
clear enhancement at small alignment angles and shape 
resonance may have been seen in the 1300-nm data, no 
evidence of shape resonance has been observed in the 
HHG spectra from 800-nm lasers. Further experiments 
dedicated to resolve this issue will be of great interest. 
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